Optimal generation of single-qubit operation from an always-on interaction 

by algebraic decoupling 
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We present a direct algebraic decoupling approach to generate arbitrary single-qubit operations 
in the presence of a constant interaction by applying local control signals. To overcome the difficulty 
of undesirable entanglement generated by the untunable interaction, we derive local control fields 
that are designed to both drive the qubit systems back to unentangled states at the end of the time 
interval over which the desired single-qubit operation is completed. This approach is seen to be 
particularly relevant for the physical implementation of solid-state quantum computation and for 
the design of low-power pulses in NMR. 



Most schemes for implementation of quantum compu- 
tation require achieving both single- and two-qubit oper- 
ations in order to realize the speed-up associated with 
quantum algorithms []]. Typically, single-qubit oper- 
ations are implemented by external local control fields 
applied to individual qubits, whereas two-qubit opera- 
tions are generated by interplay of the couplings between 
qubits and local fields [2| ■ To simplify engineering design 
and reduce decoherence channels, many physical propos- 
als of quantum computation have "always-on" and fixed 
couplings between qubits. This is often the case in ar- 
rays of solid-state qubits, which constitute a very attrac- 
tive research direction because of the inherent scalability 
of the required microfabrication techniques. For exam- 
ple, in superconducting systems the interaction between 
qubits is often coupled by a capacitor or inductor, whose 
value is fixed by the fabrication and generally cannot be 
tuned during the computations 0, y|. Several variable 
coupling schemes have been suggested 0, 0, 0, H, Ej , but 
it is generally agreed that none of these is completely sat- 
isfactory. In contrast, proposals for using electron spins 
in quantum dots do in principle allow for electrical gating 
of the exchange spin-spin interaction [Tol ITl| , but engi- 
neering such control in practice still remains extremely 
challenging (See, e.g., Ref. 01 an d references therein.) 
Other spin-coupled implementations of current interest 
include endohcdral fullcrcncs for which non-tunable mag- 
netic dipolar coupling between electron spins of neighbor- 
ing endohedral species provide the required qubit cou- 
pling mm. 

While coupling between qubits is essential for the im- 
plementation of two-qubit operations, an always-on and 
untunable coupling leads to additional complications in 
implementation of single-qubit operations, because the 
qubits may become entangled. This general issue is also 
encountered in nuclear magnetic resonance. For the short 
and high power "hard" pulses, sophisticated refocusing 
schemes have been developed to decouple the subsys- 
tems |l5j . However, for the low-power pulses often used 
in homonuclear spin systems, the radio frequency signals 
can be of the same order of magnitude as the coupling 
strengths [Til \vfy. In this "soft" pulse situation, the 



single-qubit rotations cannot be assumed to be imple- 
mented instantaneously as is the case with hard pulses, 
and consequently the interactions may strongly affect the 
intended operations. To overcome this general problem in 
the context of quantum information processing, Ref. |Tsj 
employed an encoding of logical qubits, while logic op- 
erations are performed using an analogue of the NMR 
selective recoupling method. Ref. [Tj| presented another 
encoding scheme to realize universal quantum compu- 
tation on carefully designed interaction free subspaces. 
Recently, Ref. [Ij showed that one can also avoid this 
undesirable entanglement by tuning the transition en- 
ergies of individual qubits. More complex schemes in- 
volving auxiliary degrees of freedom have also been de- 
veloped IT3 | . Here we present a direct approach for 
the design of local control signals in the presence of an 
always-on interaction. To eliminate the accompanying 
entanglement generated by this untunable interaction, 
we derive local control fields that drive the qubit sys- 
tems back to unentangled states at the end of the time 
interval over which the desired single-qubit operation is 
completed. 

We consider a two-qubit system with the following 
Hamiltonian: 
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where <j x , a y , and a z are the Pauli matrices, J the 
always-on and untunable coupling strength, u)j and <j>j 
the amplitude and phase of the external control fields, 
respectively. This Hamiltonian describes a wide range 
of two-qubit systems, e.g., solid-state quantum computa- 
tion with superconducting circuits as proposed in as 
well as NMR double-resonance J cross polarization with 
two RF fields in the XY plane of the doubly rotating 
frame |21|. Note that the Hamiltonian in Eq. JIJ can be 
transformed by local unitary operations into a Hamilto- 
nian with YY coupling and a x , <r z local terms, or into 
one with XX coupling and <7 y , a z local terms. Therefore, 
the results in this paper are also applicable to these two 
types of Hamiltonians. The spin coupling J here can be 
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neglected only for short and high-power pulses, since for 
low-power pulses u>\ and uj 2 can be of the same order 
as J. We now show how to design amplitude-modulated 
pulses that can generate an arbitrary single-qubit opera- 
tion for such a system. The control pulses are derived by 
algebraically decoupling the two-qubit Hamiltonian into 
two unentangled single-qubit systems. This approach is 
easy to implement in the physical systems of interest and 
avoids the encoding overheads associated with the encod- 
ing schemes of Refs. [T^.ll9) . 

An arbitrary local unitary operation k\ on two qubits, 
the target operation here, can be written using Eulcr's 
XYX decomposition as: 

k x = ( e - iaia * /2 (g) e- <a a°W2)( e -i/3i°y/2 ge"' ftl » /2 ) 

x(e-* n "'/ J «e- il! " r "' J ), (2) 

where ay, (3j, and 7^ are the Euler angles. In or- 
der to generate ki, we need to generate arbitrary <j x 
and a y rotations on each qubit. However, we can sim- 
plify this problem to the generation of only a x rotations 
k 2 = e~^ ia "/ 2 ® e- l72<T */ 2 from the following Hamilto- 
nian: 

Hi = ^rol + + J -<j\al, (3) 



where H\ is obtained from Eq. by setting the phases 
4>i = 02 = for external control pulses. This leads im- 
mediately to the implementations of the first and third 
operations in Eq. J2J. For the second operation in Eq. J2J), 
we can set <j>\ = <fr 2 — n/2 in Eq. JIJ to obtain a Hamil- 
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where V = e i7! / 2a */ 2 £ 
quantum operation e" 
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ate the second operation in Eq. J2J), since 



Our task has therefore been simplified to the genera- 
tion of an arbitrary local unitary operation e - l 7i<>W2 ^ 
e -272oW2 f rom ^ ne Hamiltonian Hi in Eq. ©. We now 
observe that the terms icl/2, i<J 2 /2, and ia\a 2 /2 ap- 
pearing in TJi generate the following Lie algebra: 



£ i = ol^' a x> a z a v a v a z> a v a v a z (T zi- 



(5) 



It is straightforward to show that 61 satisfies the same 
commutation relations as so(4), where so(4) denotes the 
Lie algebra formed by all the 4x4 real skew symmetric 
matrices. Therefore, ti is isomorphic to so(4). We also 
know that so(4) is isomorphic to su(2) ® su(2) To 
realize this fact, let 
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and use £ 2 to denote the Lie algebra generated by {iel, 
t€y, iel, * e x> i e y> * e *}- We have the following commuta- 
tion relations for £2: 
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It is clear that the Lie algebra t 2 satisfies the same com- 
mutation relations as i/2{a x , a y , cr z , a 2 , a 2 , a 2 }, and 
therefore it is isomorphic to su(2) ® su(2). This iso- 
morphism allows great simplification for the generation 
of single-qubit operation from Hamiltonian JSJ) , because 
it provides an algebraic way to decouple the entangling 
Hamiltonian into two unentangled single-qubit Hamilto- 
nians. To our knowledge, this fact has not been recog- 
nized, although a transformation similar to 10 was pre- 
sented in [2l|. 

We can now rewrite Hi as 

Hi = (wi - to 2 y x + J4 + (wi + uj 2 )e 2 x - Je 2 y , (7) 

and the desired generator k 2 as 

k 2 = e-^i^/a 



-1720-^/2 _ e -i((7l-72)ei + (7l+72)el) 



(8) 



We now transform the parameters in k 2 and Hi to refor- 
mulate the problem as generation of the local unitary 



-i((7i-72)^/2+(7i+7 2 )<T^/2) 
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To obtain Eqs. @ and (|10|) we simply replaced e 3 a in 
Eqs. JTJl and (JSJ) by <r 3 a /2, which is warranted by the 
fact that {iel,ie y ,iel,ie 2 ,ie 2 ,ielj,} is isomorphic to the 
Lie algebra su(2) ® su(2). Taking advantage of the com- 
mutation of a\ and cr^, the problem then naturally de- 
composes into two well-defined problems of generating 
single-qubit rotations: 

(1) Generate e - i ^~y^ll' 1 
(wi - u % )o\l2 + Jcr.i/2; and 

(2) Generate e-^+^ll 2 
{ui+u 2 )a 2 x /2~Ja 2 y /2. 

By making the transformation from £1 to su(2)(g)su(2), 
we have therefore arrived at two decoupled single-qubit 
quantum systems, solutions of which both reduce to a 
general steering problem on the Lie group SU(2) with 
dynamics determined by the Schrodinger equation: 



from the Hamiltonian 



from the Hamiltonian 



iU = (^-* x + ^ v )U, 17(0) 



(11) 
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In general, uj(t) is a time dependent external control field. 
The minimum energy control on the Lie group SU(2) has 
been studied in Ref. [2]|. Here we will use Lie-Poisson 
reduction to derive both the minimum energy and time 
optimal control. Wc will also give a good approximate 
solution to this problem. 

Pontryagin's Maximal Principle provides an important 
mathematical tool for solving the optimal steering prob- 
lem on Lie groups [24| . We solve for a control field to that 
drives the system from the initial operation C/(0) = I to a 
prescribed target operation U(T) = e~ l ' y ' 2c ' x , which also 
minimizes the cost function: 



J 



L(u)(t))dt, 



(12) 



where L(uj(t)) is a general functional of the control field 
Lo(t) and is often referred as the running cost. The first 
step is to construct the control Hamiltonian: 

H(t) = L(u(t)) + (M, -iU(tf(^a x + ia y )U(t))(13) 

where M is a constant matrix in su(2), and (X,Y) = 
Tr(XY^) is an inner product defined on su(2). For the 
ease of notation, we will suppress the time parameter 
unless otherwise specified. Proceeding further, we can 
write Eq. I|13|) as 

H = L{uj) - ]-{M,iU ] a x U)u> - ]-{M,itf<j y U)J. (14) 

We then let Pl = (M , if/1 'a x U) /2, p 2 = (M,iU i <ayU)/2, 
and p 3 = (M,iWa z U)/2. Taking the derivative of pi, 
we obtain 

px = ^(M,itfa x U) + ^(M,iU f a x U) = Jp 3 . (15) 

By similar means, we can get the dynamics for p 2 and 
p 3 . Combining these together, we have the dynamics for 
Pj as 



Pi = JP3 
p 2 = ~UP3, 

i>3 = ^Pi - Jpi- 



(16) 



Eq. (|16H is indeed the Lie-Poisson equation on SU(2) |25| . 
It is easy to verify that the functions C\ = p\ + p\ + p\ 
and C 2 = Pi + 2Jp 2 are both invariant along the sys- 
tem trajectory. For a general optimal control problem 
on Lie group SU(2), one usually needs to solve a set of 
six differential equations. However, here we can obtain 
the reduced dynamics of three differential equations in 
Eq. (|Pd|| by using the Lie-Poisson reduction theorem (25[ . 
The two conditions for this theorem are the vector field is 
right invariant and the cost function being independent 
of the quantum states, and they are both satisfied in this 
case. This reduces the number of the differential equa- 
tions to be solved by half. This approach is a general 



technique applicable also to analysis of optimal control 
on multi-qubit systems. 

Now we derive the time optimal control, that is, the 
form of control field lu that achieves the target operation 
in the minimum time possible. In this case, the cost 
function can be written as 



min J 



Idt, 



and the running cost L(lu) = 1 is independent of control 
field lu. The corresponding control Hamiltonian is 



H = 1 - LOpi - Jp2 



(17) 



Let us assume that the control field to is restricted to 
an interval [u> m , u>m) due to physical bounds on the avail- 
able frequencies. From Pontryagin's Maximum Principle, 
the optimal control field ui{t) will minimize the control 
Hamiltonian (|17|) pointwisc. Therefore, u)(t) can take 
only extremal values: 



uj m , \ip\{t) > 0, 
ojm, Hpi(t) < 0. 



This means that the time optimal control will switch back 
and forth between two extremal control values. This pro- 
cedure is usually called Bang-Bang control in the control 
theory literature. (Note that the term "Bang-Bang con- 
trol" has recently been adopted with a somewhat differ- 
ent meaning in the study of dynamical coupling of open 
quantum systems [2(| |23, , where it is referred to 

performing a set of instantaneously or as fast as physi- 
cally possible hard pulses.) In Ref. @, we have derived 
a constructive approach to achieve a desired single-qubit 
target quantum operation by switching between two such 
constant control fields. From the arguments above, we 
conclude now that this is indeed a time optimal control 
strategy. 

Next we consider the minimum energy control which 
is encoded by the following cost function: 



min J 



1 



\t)dt. 



The control Hamiltonian is now 
1 



n 



-U! - UJpi - Jp 2 



and thus the minimum energy control is Lo(t) = pi(t). 
Taking the derivative of p\ in Eq. IjRjfl . we have 



[ T J )~ k'- 



(18) 



where Ci = p\ + 2Jp 2 is constant along the optimal tra- 
jectory. The solution of this differential equation is given 
by the Jacobi elliptic function Cn [23j : 



Q(t) = 2bkCn(bt + f,k), 



(19) 
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where b, /, and k are real numbers. We can numerically 
determine the values of these parameters such that the 
control function l(19|) steers the Hamiltonian from the 
initial operator U(0) — I to the target operation U(T) = 

e -i7/2(T x 

Finally we give an approximate solution to the con- 
trol problem in terms of a sinusoidal control function 
u> = Acos(vt), and show that this can be further fidelity 
optimized. Letting U\ = e m / 2ayt U , we transform the 
system <|11[1 into the following form: 

if/i = {4(1 + cos2u£)y + (J - v)y + ^sin2wiy}C/ 1 

(20) 

The usual technique to solve this differential equation in 
the context of NMR is to drop the oscillating terms to get 
a time-independent Hamiltonian (rotating wave approx- 
imation). In order to get a more accurate solution, we 
use the Wei-Norman formula to transform Eq. I(20|) 
to a dynamical system on R 3 . From Euler's ZXZ decom- 
position, a general solution to Eq. (|20|l can be written 
as 

U\(t) = e - iQ i( t )^/2 e -ia2(t)o' :c /2 e -iQ3(t)o' i /2^ ^1) 

Taking the time derivative of U\, we have 

iUi = |(d2 cosai + d3 sinai sin a -2)-^ + («2 sinai 
— 6(3 cosai sin 0:2)-^ + az cosc^)-^} Ui- (22) 

Comparing Eqs. lf2*U|) and (J55Jl, we get 



A/2(l + cos2vt) 
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Therefore, we obtain a set of differential equations of the 
parameters ctj in M. 3 : 







sin Qi 
tan Q2 


COS Qi 

tan oli 






cos ai 
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d 3 




sin Qi 
_ sin 02 


COS OL\ 

sin a.2 



A/2(l + cos2vt) 
J -v 
A/2 sin 2vt 



.(23) 



Rigorous solution of Eq. 1)23(1 cannot be given in terms of 
elementary functions. However, when v is close to J, a 
very good approximate solution can be found as: 



ai(t) 

a 2 (t) 
a 3 (t) 



A . 2 
— sin vt, 
2v 



= 0. 



Hence a sinusoidal control function u> — Acos(vt) can 
achieve the following unitary operation at a final time T: 



To generate an operation e~ ll / 2(7x with 7 S [0, 27r], we 
only need that v = J, AT/ 2 = 7, and JT = 2nir, where 
n is an integer. This leads to the conditions: 



lJ „ 2nn 
v = J, A=^, T=—. 

nir J 



(25) 



From Eq. JHJ, to implement a two-qubit local unitary 
operation e~ nia!C ^ 2 ® e _l72CTx / 2 , we require that 



COS Jt, Ul-2 



I2J 



COS Jt, 



(26) 



with a time duration 2nn/J. Note that this is an ap- 
proximate solution. Numerical simulations reveal that 
the greater the value of n, the better the approximation. 
To improve the accuracy, we can numerically search the 
variable parameters that maximize the fidelity of the ac- 
tual achieved operation and the desired target operation, 
using parameter values given by Eq. I|25|) as a starting 
guess. This gives rise to a fidelity optimized control. 

We illustrate the control strategies aforementioned 
with an example relevant to NMR experiments on 
homonuclear spin systems |0, ^} . Let J = 200 Hz in 
Eq. Q), and our target the generation of a 90° rotation 
about x-axis on the first qubit, that is, the operation 
Choosing (j>x = fa = 0, n = 1, the analytic 
approximation Eq. 1251) yields 

v = 200 Hz, A = 100 Hz, T = 10tt ms, (27) 

which leads to an approximate solution: 

u)i = 100cos200i, lo 2 = 0. (28) 

Numerical optimization via the maximization of the fi- 
delity of actual achieved operation and the desired target 
operation, using the parameters in Eq. I|27|) as an initial 
guess, leads to the following fidelity optimized control 
function 



wx = 98.062 cos 196.900t, oj 2 = 0, 



(29) 



U(T) 



e tvT* y /2 e -iai (T)<r 2 /2 e -ia 2 (7> x /2 ^A) 



with a final time T — 31.911 ms. These two control 
functions, the analytic approximation and the fidelity op- 
timized function, are shown in Fig.^A). The alternative 
minimum energy control function is determined numeri- 
cally, resulting in parameters b = 210.744, / = —0.00549, 
and k = 0.00236, in Eq. i|19fl . and the minimum energy 
control: 

ui = 99.678Cn(210.744t - 0.00549, 0.00236). (30) 

This minimum energy control solution is seen to be very 
close to the fidelity optimized control function of Eq. 12911 . 
The difference between these two control functions is 
plotted in Fig.QlB). 

In summary, we have presented a design of direct con- 
trol pulses that implement any arbitrary single spin op- 
eration in a quantum system with an always-on interac- 
tion. This is crucial for the physical implementation of 
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FIG. 1: Control functions that generate the single spin oper- 
ation e~ l7r/4rT x when J = 200 Hz. (A) Dashed line: approxi- 
mate control in Eq. 1281 : solid line: fidelity optimized control 
in Eq. H29II : (B) the difference between fidelity optimized con- 
trol and minimum energy control in Eq. 13011 . 

solid-state quantum computations and for the low-power 
pulses design in NMR. We used an algebraic approach 



to decouple the two-qubit system to overcome the diffi- 
culty of undesirable entanglement generated by the un- 
tunablc interaction. To generate the desired target quan- 
tum operations, we derived three control strategies, i.e., 
minimum time, minimum energy, and an approximate 
solution. The first two optimal strategies implement the 
desired target quantum operation exactly, whereas the 
approximate control can be further optimized for perfect 
fidelity. The advantage of the approximate control is that 
we can easily determine the control parameters, and it is 
close to the minimum energy control. 
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